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Abstract 

Using numerical simulations, we observe phase locking, Arnold tongues, and 
Devil's staircases for vortex lattices driven at varying angles with respect to 
an underlying superconducting periodic pinning array. This rich structure 
should be observable in transport measurements. The transverse V{I) curves 
have a Devil's staircase structure, with plateaus occurring near the driving 
angles along symmetry directions of the pinning array. Each of the plateaus 
corresponds to a different dynamical phase with a distinctive vortex structure 
and flow pattern. 
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Introduction. — Numerous nonlinear driven systems in physics, astronomy, and engineer- 
ing exhibit striking responses with complex phase-locking plateaus characterized by Devil's 
staircases, Arnold tongues, and Farey trees Here, we present the first evidence that 

these structures can be observed in bulk superconductors. 

Driven vortex lattices (VLs) interacting with either random or periodic disorder have 
attracted growing interest due to the rich variety of nonequilibrium dynamic phases which 
are observed in these systems. These phases include the elastic and plastic fiow of vortices 
which can be related to VL order and transport properties Periodic pinning arrays 

interacting with VLs are now attracting increasing attention as recent experiments with 
patterns of holes and magnetic dots |Tl| have produced interesting commensurability 
effects and enhanced pinning. These systems are an excellent realization of an elastic lattice 
interacting with a periodic substrate that is found in a wide variety of condensed matter 
systems including charge- density- waves, Josephson-j unction arrays, and Frenkel-Kontorova- 
type models of friction (see, e.g., ll^l)- An interesting aspect of periodic pinning arrays 
that has not been addressed so far is how the symmetry properties of the array affect the 
transport properties as the VL is driven at different angles. 

We find that as a slowly increasing transverse force is applied to a VL already moving in 
the longitudinal direction, the VL undergoes a remarkable series of locking transitions that 
significantly affect both the VL ordering and transport properties. These locking phases 
occur when the direction of the vortex motion locks with a symmetry direction of the pin- 
ning array. As the VL passes through these phases, the transverse velocity component as 
a function of increasing transverse drive shows a series of plateaus which form a Devil's 
staircase structure At the boundaries of certain locked phases the VL undergoes a 

transition to a plastic flow phase in which defects are generated in the VL. In the locked 
phases the VL undergoes elastic flow in static ID channels and the overall VL has a variety 
of orderings, including triangular and square. 

Simulation. — We consider a 2D slice of A'^ 3D rigid vortices interacting with a square 
array of Np parabolic wells, with lattice constant a, and periodic boundary conditions. We 
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integrate the equations of vortex motion fj = f^^^ + + = rjVi. The total force fj 

on vortex i includes interactions with other vortices , pinning f^^ by parabolic wells, and 

an applied driving force = /^^x. The vortex-vortex interaction between vortex i and the 

other vortices is f™ = Y^fZi foKi{\ri — rj|/A)fij, where Ki{r/X) is a modified Bessel 

function, A is the penetration depth, /o = ^q/Stt^A^, = (rj — rj)/|rj — rj|, and we set 

Tj = 1. All lengths, fields, and forces are given in units of A, $0/-^^; and /o, respectively. For 

most of the results presented here the number of vortices is close to the number of pinning 

sites, = 1.062A''p. We have conducted a series of simulations with different pinning 

parameters so that accurate phase diagrams of the dynamic phases can be obtained. In order 

to investigate finite size effects we have examined system sizes varying between 36A x 36A 

and 108A x 108A, with A^^ between A^^ = 550 and A'„ = 4955. 

Voltage-Current Response. — First, the VL ground state at zero applied driving force 

is found by simulated annealing (i.e., by cooling the VL from high-T). After a low energy 

ground state is found, a slowly increasing driving force, f^, is applied along the horizontal 

symmetry axis of the square pinning. We find that increasing /j, in increments of O.OOl/o 

every 400 MD steps, from /x = to = 3.0/o, is slow enough that the vortex dynamics 

does not depend on the rate of increase of f^^. Once is brought to 3.0/o it is held constant 

while a force, which we label fy, is applied in the transverse or y-direction. We increase 

fy from to 3.25/0, also in increments of O.OOl/o every 400 MD steps. The total driving 

force has a net magnitude of fd = {fl + fyY^"^ at an angle 9 = tan^^ (fy / f^) with respect 

to the x-direction. We compute the average velocity of the moving vortices in both the 

longitudinal = (l/A"!,) Yli^i^i " ^ and the transverse Vy = (l/Ny) Ylj^=i^i " y direction, 

as fy is increased. Velocity versus driving plots correspond to experimentally measurable 

volt age- current V{I) curves. 

In Fig. 1(a) we present a typical plot of and Vy. For fy < 0.4/o, Vy = indicating 

that the VL is pinned in the y-direction even though the VL is moving in the x-direction. 

Depinning in the transverse direction occurs at fy = 0.4/o, as indicated by the sharp jump up 

in Vy. We label this critical transverse depinning force fy. A jump up in Vx is also observed 
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at fy. As fy is linearly increased, Vy does not grow linearly but instead in a remarkable series 
of jumps and plateaus of varying sizes 0] . Along the plateaus Vy is constant or increasing 
very slowly, indicating that the vortex motion is locked in a certain direction for a finite range 
of increasing fy. The small jumps and dips in Vx correspond to the onset of plateaus in Vy. 
The plateaus in Vy occur when the ratio of fy to f^ is near a rational value: fy/ fx = p/q, 
where p and q are integers. In Fig. 1(a) the largest plateaus occur at p/q = 0, 1/3, 1/2, 2/3 
and 1. Fig. 1(b) shows a blow-up of a region in Fig. 1(a) for values of fy = O.6/0 to 2.1/o, 
where additional plateaus at p/q = 1/5,1/4,2/5,3/7, and 3/5 are highlighted. For larger 
system sizes we find exactly the same behavior in Vy and Vx as observed in Fig. 1, indicating 
that it is independent of the system size. 

Vortex Dynamics and the Origin of the Plateaus. — To understand why the plateaus 
occur as well as the VL dynamics in the plateau and non-plateau regions, in Fig. 2(a-d) we 
plot the vortex trajectories for rational ratios of fy/fx = 0, 1/2, 1, and the irrational ratio 
fy/ fx = 27r/ll = 0.571... In Fig. 2(a), where fy < fy, the vortex motion traverses pin sites 
periodically and it is only along the x-direction — with the vortex flow restricted in ID paths 
along the pinning rows. This periodic ID motion persists up to fy = fy, at which point 
the vortices also begin to flow in the y-direction. In Fig. 2(b), for p/q = 1/2 where a large 
plateau in Vy is observed in Fig. 1, the vortices again exhibit periodic motion and flow in ID 
channels along the pinning sites — and along a symmetry axis of the pinning array at an angle 
6 = tan~^(l/2) from the x-axis. A similar periodic ID motion is seen in (d) for fy/fx = 1, 
with the VL motion at 45° from the x-axis. In Fig. 2(c), at the irrational fy/fx ratio, the 
vortex trajectories are different than those observed in Fig. 2(a,b,d). Here the quasiperiodic 
vortex trajectories drift over time, eventually covering the sample (i.e., ergodic-like motion). 
In general, the plateau regions (with rational fy/fx) in K/(-^) correspond to periodic ID 
vortex trajectories, while the non-plateau regions produce quasiperiodic trajectories [p],|2|,p!3| . 

To understand how the vortex motion locks into certain driving angles, we first consider 
the case fy/fx = 0. Here the vortices move along the pinning rows in ID paths, with 
each vortex traversing a distance a — 2rp between pinning sites, as seen in Fig. 2(a). An 
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application of a transverse force fy causes the moving vortices to drift a small distance in the 
y-direction. Once the vortices interact with the pinning sites, they feel a force that moves 
them towards the center of the pinning site which keeps them locked along the x-direction. 
When fy is large enough, fy ^ fx tan(rp/a) |jl3[, the vortices are able to break off from 
moving only along the x-direction and start moving in the ?/-direction as well. 

As fy is increased beyond fy, the net driving force vector will be at an angle with 
the horizontal. Due to the symmetry of the square pinning array, along the angles where 
9 = tan^^ (p/q), the vortices encounter pinning sites periodically spaced a distance ae apart. 
This distance is related to the pinning lattice constant a by = a(j9^ + g^)^/^. Along these 
commensurate angles, the vortex motion will be periodic and locked in ID channels in a 
similar manner as the fy/fx = case. The force needed to depin the vortices from the 
commensurate angles will vary since varies. For values where ag is small, the vortices will 
move only a small distance between pinning sites, so a higher depining force is needed. For 
large ag the vortices will move a much longer distance before encountering the pinning sites, 
so a much smaller depinning force is needed. This is in agreement with Fig. 1 where the 
largest plateaus (due to enhanced pinning) occur for values of p/q that produce the lowest 
distance between pinning sites, that is the smallest ag (i.e., p/q = 0/1, 1/1, and 1/2). 

The onset of certain plateaus coincide with a variety of structural transitions in the 
VL. We quantify this angle-dependent evolving topological order by using the Voronoi (or 
Wigner-Seitz) construction to obtain the fraction of vortices with coordination numbers six, 
Pq, and four, P4. In Fig. 3(a,b) we show the evolution of Pq and P4 as fy is increased, for 
the same system as in Fig. 1. For fy < fy, Pq ^ 0.68, indicating a mostly triangular VL. 

fy = fy) ^ dip in Pq, along with direct observation of the VL flow, show that the VL 
disorders due to plastic deformations. Right after the initial dip in Pg the VL suddenly 
regains considerable triangular ordering, as indicated by Pg ~ 0.95. Small dips in Pg can 
be seen near the 1/4, 1/3, and 2/3 locking regions. At the 1/2 locking region the VL 
is considerably disordered, as indicated by the sharp drop in Pg. This is consistent with 
Fig. 3(c), where both the vortex positions and Voronoi polygons are shown for a 12A x 12A 
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region in the 1/2 locking region. At the 1/1 locking region Pq drops almost to zero while 

P4 increases to about 0.9, indicating a structural phase transition from a triangular to a 

square VL. Here, the fy = fx symmetric drive is what produces a moving square VL. The 

less symmetric drives {2fy = and 3fy = fx), produce more distorted squares |T3|. For the 

special case when fy = and for the B used in Fig. 3, correlations between nearby VL rows 

are strong, and near 2/3 of the VL has triangular order (which diminishes for weaker 5's). 

In Fig. 3(d,e) the vortex positions and Voronoi polygons are shown for (d) right before the 

transition to the 1/1 locking region and (e) in the 1/1 locking region showing the triangular 

and square ordering of the VL respectively. Right at the boundaries of the 1/1 phase, the 

VL is strongly disordered and has a similar structure to Fig. 3(c). 

Phase Diagrams with Arnold Tongues. — We have derived five phase diagrams which 

indicate the evolution of the plateau regions versus the following parameters: fp, Up, r^, 
commensurability, and disorder. These five phase diagrams are all very similar, and thus 
here we present only one: Fig. 4(a). This is obtained by conducting a series of simulations 
in which the maximum pinning force fp is varied between 0.25 < fp/fo < 2.75. The 
phase diagram shows 18 clearly defined (shaded) Arnold tongues or plateaus [l|-§]. As fp 
is decreased the widths of the tongues also show a corresponding decrease. For fp/fo > 2.5 
several locking phases are lost (i.e., 1/6, 4/7, 5/6) due to overlapping by other locking 
regions. For fp/fo < 1.0 only the strongest plateau regions can be resolved within the 
accuracy of our calculations. 

The phase diagram in Fig. 4(a) has the same structure as Arnold tongues |]1]-§1 found 
in phase locking systems where the widths of the tongues, or locking regions, increase as 
the nonlinear coupling increases. Here, the coupling is between the vortices and the pinning 
array, and is increased with increasing fp, Tp, Up, density of vortices (i.e., the commensu- 
rability B/B^), and pin-location order [|T^. In Fig. 4(b) we present the width of the 0/1 
locking region for varying pinning density in units of the pinning lattice constant a. As a 
decreases the width of the locked region increases. This can be understood by considering 
that as a decreases the vortices in the locked region will move a smaller distance between 
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pinning sites; thus a higher transverse force is needed to break the vortices away from the 
locked region. The widths of the other locked regions show the same behavior as the 0/1 



region for increasing a [13 



We have also examined the effects of pin disorder on the width of the locking regions by 
conducting a series of simulations in which the pinning sites are randomly displaced up to 
an amount 5r away from the perfectly square pinning lattice. We consider the case where 
(5r = a/2 to be a good approximation to a random pinning array. In Fig. 4(c), we examine 
how the width of the 0/1 locking region, fy, decreases as 5r is increased. It is of interest 
to compare our results for large disorder with Ref [§](a) in which a nonzero transverse 
critical force fy was predicted. Recent T = MD simulations have observed extremely 
small transverse barriers 0(b). We find that for large disorder, 5r = a/2, a true transverse 
barrier (i.e., = 0) is not observed. Also, for a triangular array of pins, the plateaus occur 
for e = tan~i(V3p/(2g + 1)). 

Summary. — In conclusion, we have found that as an increasing transverse force is applied 
to a strongly driven VL interacting with a periodic pinning array, the VL undergoes a 
remarkable series of locking transitions in which both the VL order and flow patterns change. 
As the VL passes through these transitions, Vy exhibits a striking series of plateaus forming 
a Devil's staircase structure. The width variations of these plateaus with different pinning 
form Arnold tongues which can be indexed via a Farey tree construction. These locking 
effects occur whenever the VL is driven along a symmetry angle of the pinning array. For a 
square pinning array, the locking phases occur when driving in the longitudinal direction is 
a rational ratio, fy/fx = p/q- These predictions can be tested experimentally and we hope 
that this work will motivate several novel experiments. Moreover, other candidate systems 
where these predictions may be accessible include: driven Wigner crystals interacting with 
a periodic array of donors, driven colloids interacting with optical-trap arrays, spin- and 
charge- density waves, Josephson-j unction arrays, and solid friction experiments. 

We thank F. Marchesoni, M. Bretz, and very specially C.J. Olson for useful discussions. 
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FIGURES 



FIG. 1. (a) Average longitudinal Vx (upper curve) and transverse Vy velocities versus the 
transverse driving force fy for a 36A x 36A sample with a square pinning array, desity of field lines 
B satisfying B/B^ = 1.062, matching field B^ = 0.4$o/A^, density of pinning sites Up = 0.4/A^, 
fp = 2.5/0, a = 1.57A, and Vp = 0.3A. fx is fixed at fx = 3.0/o. Plateaus are seen in Vy near values 
where fy/ fx = p/q, where p and q are integers. The largest plateaus (at 0/1, 1/3, 1/2, 2/3, and 
1/1) are clearly seen, (b) shows a blow up of Vy from (a), for fy = O.6/0 to 2.1/o, where additional 
plateaus at 1/5, 1/4, 2/5, 3/7, and 3/5 can be seen more clearly. The overall structure in Vy is 
that of a Devil's staircase [1-3] . 



FIG. 2. The vortex trajectories for a subset of the system in Fig. 1 at the plateau re- 
gions (a) fy/fx = 0, (b) fy/fx = 1/2, and (d) fy/fx = 1/1, and at the non-plateau region (c) 
fy/fx = 27r/ll = 0.571... At the plateau regions, the vortices move in ID channels, periodi- 
cally along the pinning rows; while at the non-plateau regions the vortices exhibit quasiperiodic 
trajectories. 
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FIG. 3. The fraction of (a) six-fold Pq and (b) four-fold P4 coordinated vortices versus trans- 
verse driving force fy, for the same system as in Fig. 1. Large drops in Pq can be seen at fy, as 
well as at the 1/2 and 1/1 locking regions. Smaller dips in Pq can be seen at the 1/4, 1/3, and 2/3 
plateaus regions. At the 1/1 transition P4 rises to 0.9 indicating a transition to a square VL. In 
(c,d,e) both the vortex positions and Voronoi polygons for a subset of the VL can be seen for: (c) 
the 1/2 locking region, where a disordered VL is observed; (d) right before the 1/1 plateau, with 
a triangular VL; and (e) at the 1/1 plateau, with a square VL. 

FIG. 4. (a) Phase diagram, for the system in Fig. 1, showing Arnold tongues (shaded); i.e., 
the widths of the locking regions versus fp. As fp decreases, the tongues (locking regions with 
periodic trajectories) shrink. In (b,c) the width of the first locking region or fy is shown versus 
pinning lattice constant a (b) and disorder 6r/2a (c). 
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